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Abstract 

The multicommodity flow problem is NP-hard already for two commodities 
over bipartite graphs. Nonetheless, using our recent theory of n-fold integer 
programming and extensions developed herein, we are able to establish the 
surprising polynomial time solvability of the problem in two broad situations. 

1 Introduction 

The multicommodity transshipment problem is a very general flow problem which 
seeks minimum cost routing of several discrete commodities over a digraph subject to 
vertex demand and edge capacity constraints. The data for the problem is as follows 
(see Figure [l] below for a small example). There is a digraph G with s vertices and 
t edges. There are / types of commodities. Each commodity has a demand vector 
d!' G with (ij the demand for commodity k at vertex v (interpreted as supply 
when positive and consumption when negative). Each edge e has a capacity 
(upper bound on the combined flow of all commodities on it). A multicommodity 
transshipment is a vector x = {x^, . . . , x') with x'' G Z*,^ for all k and the flow of 
commodity k on edge e, satisfying the capacity constraint Ylk=i — ""e for each 
edge e and demand constraint X]ee5+(j;) -^e ~ J2ees-{v) = each vertex v and 

commodity k (with 6^{v), S~{v) the sets of edges entering and leaving vertex v). 

The cost of transshipment x is defined as follows. There are cost functions 
/ejS'g : Z ^ Z for each edge and each edge-commodity pair. The transshipment 
cost on edge e is feiYlk=i^e) + Ylk=i 9ei^e) with the first term being the value of 
/e on the combined flow of all commodities on e and the second term being the sum 
of costs that depend on both the edge and the commodity. The total cost is 

e=l V \k=l / k=l / 
* Supported in part by a grant from ISF - the Israel Science Foundation 
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Multicommodity Transshipment Example 

Data: 

digraph G 

two commodities: red and green 
edge capacities Ug ulimited 

edge costs ^x\+x\)\={x\+x%f and Q\{x\y.=g%{x%y.=0 
vertex demands: 

di := (3 -1 -2) 

d2:=(-3 2 1) 

Solution: 

Xi = (3 2 0) 
X2= (0 2 3) 

Cost: 

(3+0)2+(2+2)2+(0+3)2=34 




Figure 1: Multicommodity Transshipment Example 



Our results apply to cost functions which can be standard linear or convex such 
as ae\ Ylk=i ^el^" + Ylk=i 7e kel^^ some nonnegative integers ae, Pe,7e^ S^, which 
take into account the increase in cost due to channel congestion when subject to 
heavy traffic or communication load (with the linear case obtained by jSe = (5g=l). 

The problem is generally hard: even deciding if a feasible transshipment exists 
(regardless of its cost) is NP-complete already in the following two very special 
cases: first, with only / = 2 commodities over the complete bipartite digraphs Km,n 
(oriented from one side to the other) [H |5]; and second, with variable number of 
commodities over the digraphs K^ n with m = 3 vertices on one side (see Section 4). 

Nonetheless, using the theory of n-fold integer programming recently introduced 
in [21 El E] and extensions developed herein, we are able to establish the surprising 
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polynomial time solvability of the problem, with either standard linear costs or more 
general costs with nonlinear convex functions /e, g^, in two situations as follows. 

First, over any fixed digraph, we can solve the problem with a variable number / 
of commodities (hence termed the many- commodity transshipment problem). This 
problem may seem at a first glance very restricted: however, even for the single tiny 
bipartite digraph K^ ^, we are not aware as of yet of any solution method other than 
the one provided herein; and as noted, the problem is NP-hard for the digraphs 
K^^n- Our first theorem is the following (see Section 3 for the precise statement). 

Theorem 1.1 For any fixed digraph G, the (convex) many- commodity transship- 
ment problem with variable I commodities over G can be solved in polynomial time. 



We also point out the following immediate corollary of Theorem 1.1 



Corollary 1.2 For any fixed s, the (convex) many- commodity transshipment prob- 
lem with variable I commodities on any s-vertex digraph is polynomial time solvable. 



The complexity of the algorithm of Theorem 1.1 involves a term of O (l^^^^^ where 
g{G) is the Graver complexity of G, a fascinating new digraph invariant about which 
very little is known (even g{K3^4) is as yet unknown), see discussion in Section 4. 

Second, when the number / of commodities is fixed, we can solve the problem 
over any bipartite subdigraph of Km,n (the so-called multicommodity transportation 
problem) with fixed number m of suppliers and variable number n of consumers. 
This is very natural in operations research applications where few facilities serve 
many customers. Here each commodity type k may have its own volume Vk per unit. 
Note again that if I is variable then the problem is NP-hard already for m = 3, so our 
following second theorem is best possible (see Section 3 for the precise statement). 

Theorem 1.3 For fixed I commodities and m suppliers, the (convex) multicommod- 
ity transportation problem with variable n consumers is polynomial time solvable. 

We point out that the running time of our algorithms depends naturally on the 
binary- encoding length {d'^,Ue) of the numerical part of the data consisting of the 
demands and capacities (see Section 3), so our algorithms can handle very large 
numbers. To get such polynomial running time even in the much more limited 
situation when both the digraph and the number of commodities are fixed (where 
the number It of variables becomes fixed) and where the cost functions are linear, 
one needs off-hand the algorithm of integer programming in fixed dimension [TT] . 
However, Theorems |1. 1| and 1.3 involve variable dimension and [TT] does not apply. 
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In Section 2 we review the recent theory of ra-fold integer programming and es- 
tabhsh a new theorem enabhng the solvabihty of a generahzed class of ra-fold integer 
programs. In Section 3 we use the results of Section 2 to obtain our multicommodity 
flow Theorems 11.11 and 11.31 We conclude in Section 4 with a short discussion. 



2 A^-fold integer programming 
2.1 Background 

Linear integer programming is the following fundamental optimization problem, 

min {wx : x G Z" , Ax = b , I < x < u} , 

where A is an integer m x n matrix, b G Z"^, and /, m G with Zoo := Z l±) {±oo}. 
It is generally NP-hard, but polynomial time solvable in two fundamental situations: 
the dimension is fixed ^llj ; the underlying matrix is totally unimodular [lOj . 

Recently, in [2], a new fundamental polynomial time solvable situation was dis- 
covered. We proceed to describe this class of so-termed n-fold integer programs. 

An (r, s) X t himatrix is a matrix A consisting of two blocks Ai, A2, with Ai its 
r X t submatrix consisting of the first r rows and A2 its s x t submatrix consisting 
of the last s rows. The n-fold product of A is the following (r + ns) x nt matrix. 





(A, 


Ai ■ 


■ A, 






A2 


■ 


■ 









A2 ■ 


■ 






I 


• 


• A2 


J 



The following result of [2] asserts that n-fold integer programs are efficiently solvable. 



Theorem 2.1 [2] For every fixed integer (r, s) xt himatrix A, there is an algorithm 
that, given positive integer n, w & Z"*, b G Z'+"**, and l,u E Z^, solves in time 
which is polynomial in n and in the binary- encoding length {w,b,l,u) of the rest of 
the data, the following so-termed linear n-fold integer programming problem, 

min [wx : x G Z"* , A^'^^x = 6 , I < x < u] . 

Some explanatory notes are in order. First, the dimension of an n-fold integer 
program is nt and is variable. Second, n-fold products A^"'^ are highly non totally 
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unimodular: the ra-fold product of the simple (0, 1) x 1 bimatrix with Ai empty and 
A2 := 2 satisfies A^"'^ = 2/„ and has exponential determinant 2". So this is indeed a 
class of programs which cannot be solved by methods of fixed dimension or totally 
unimodular matrices. Third, this class of programs turns out to be very natural 
and has numerous applications, the most generic being to integer optimization over 
multidimensional tables. In fact it is universal: the results of [5] imply that any 
integer program is an n-fold program over some simple bimatrix A, see Section |4j 

The above theorem extends to n-fold integer programming with nonlinear objec- 
tive functions as well. The following two results, from [5] and P] respectively, assert 
that the maximization and minimization of certain convex functions over n-fold in- 
teger programs can also be done in polynomial time. The function / is presented 
by a comparison oracle that for any two vectors x,y can check if f{x) < f{y)- 

Theorem 2.2 [3J For every fixed d and (r, s) x t integer bimatrix A, there is an 
algorithm that, given n, bounds l,u & Z^, integer d x nt matrix W, b G Z^'^^'^ , 
and convex function / : Z'^ — M presented by a comparison oracle, solves in time 
polynomial in n and {l,u, W, b) , the convex n-fold integer maximization problem 

max{f{Wx) : a; G Z"* , ^("^x = 6 , I < x < u} . 

In the next theorem, / is separable convex, namely f{x) = Ylifii^i) with each 
fi univariate convex. The running time depends also on log / with / the maximum 
value of |/(x)| over the feasible set (/ is not needed to be part of the input). 

Theorem 2.3 [9J For every fixed integer (r, s) xt bimatrix A, there is an algorithm 
that, given n, lower and upper bounds l,u & Z^, b G Z^"*""*, and separable convex 
function f : Z"* —>■ Z presented by a comparison oracle, solves in time which is 
polynomial in n and {l,u, b, f) the convex n-fold integer minimization problem 

min{/(x) : a; G Z"* , A^''^x = b, l<x<u] . 



2.2 Generalization 



We now provide a broad generalization of Theorem 2^ which will be useful for the 
multicommodity fiow applications to follow and is interesting on its own right. 

We need to review some material from |21 E]- We make use of a partial order □ 
on defined as follows. For two vectors x,y E M" we write x ^ y if xiyi > and 
l^^il ^ \yi\ for i = 1, . . . ,n, that is, x and y lie in the same orthant of MJ^ and each 
component of x is bounded by the corresponding component of y in absolute value. 
A classical lemma of Gordan |7] implies that every subset of Z" has finitely-many 
□-minimal elements. The following fundamental object was introduced in [8j. 
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Definition 2.4 The Graver basis of an integer matrix A is defined to be the finite 
set G{A) C Z" of □-minimal elements in {s G Z" : Ax = 0, x 7^ 0}. 

The Graver basis is typically exponential and cannot be written down, let alone 
computed, in polynomial time. However, we have the following lemma from [2]. 

Lemma 2.5 For every fixed integer bimatrix A there is an algorithm that, given n, 
obtains the Graver basis ^(A^"^) of the n-fold product of A in time polynomial in n. 

We also need the following lemma from P] showing the usefulness of Graver bases. 

Lemma 2.6 There is an algorithm that, given an integer mxn matrix A, its Graver 
basis Q{A), l,u & ZJ^, b G Z'", and separable convex function / : — Z presented 
by a comparison oracle, solves in time polynomial in {A, Q{A), I, u, b, f), the program 

min{/(x) : x G Z" , Ax = b , I < x < u} . 



Note that Lemmas 2^ and 2^ together imply at once Theorem mentioned above. 
We proceed with two new lemmas needed in the proof of our generalized theorem. 



Lemma 2.7 For every fixed integer (r, s) xt bimatrix A and {p, q) xt bimatrix W , 
there is an algorithm that, given any positive integer n, computes in time polynomial 
in n, the Graver basis Q{B) of the following {r + ns+p + nq) x (nt + p + nq) matrix, 



B 



Proof. Let D be the {r + p, s + q) x {t + p + q) bimatrix whose blocks are defined by 







Ip 



Do 



A2 
W2 /, 



Apply the algorithm of Lemma |2.5| and compute in polynomial time the Graver 
basis ^(D*^")) of the n-fold product of D, which is the following matrix: 



/ Ai 
Wi Ip 


Ai 
Wi Ip 




\ 

Wi Ip 


^2 

W2 Ig 
















A2 

W2 Ig 

















\ 








^2 
W2 /, / 
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Suitable row and column permutations applied to D^'^^ give the following matrix: 



C := 



Ai 


Ai ■ 


A 

■ Ai 





• 








• 





A2 


• 


■ 





■ 


■ 





■• 


■ 





A2 ■ 


• 





• 


• 





• 


• 





• 


• A2 





• 


• 





• 


• 


Wi 


Wi ■ 


■ Wi 


Ip 


Ip ■ 


• Ip 





• 


• 


W2 


• 


■ 





• 


■ 


I, 


• 


• 





W2 ■ 


• 





• 


• 





I, ■ 


• 





• 


• W2 





• 


• 





• 


• I,, 



Obtain the Graver basis G{C) in polynomial time from Q^D^"^^) by permuting the 
entries of each element of the latter by the permutation of the columns of ^(D^")) 
that is used to get C (the permutation of the rows does not affect the Graver basis) . 

Now, note that the matrix B can be obtained from C by dropping all but the first 
p columns in the second block. Consider any element in Q{C), indexed, according to 
the block structure, as (a;^, a;^, . . . , x", y^.y"^, . . . , y", z^, z"^, . . . , z^). Clearly, if = 
for A; = 2, . . . , n then the restriction (x^, x^, . . . , y^, z^, z"^, ■ ■ ■ , -2") of this element 
is in the Graver basis of B. On the other hand, if (x^, x^, . . . , x", y^, z^, z'^, . . . , z^) 
is any element in Q{B) then its extension (x^, x^, . . . , x", y^, 0, . . . , 0, . . . , 

is clearly in Q{C). So the Graver basis of B can be obtained in polynomial time by 



, X , y ,2; , 



g{B) {{x\ 
This completes the proof. □ 



: (x^...,x^y^o,...,o,^^...,^")e^;(c)}. 



In the next lemma and theorem, as before, / and g denote the maximum values 
of |/(l^x)| and |5'(x)| over the feasible set (/, g do not need to be part of the input). 

Lemma 2.8 There is an algorithm that, given an integer mxn matrix A, an integer 
dxn matrix W, l,u e Z^, 1, u e Z^, b e Z"*, the Graver basis Q{B) of 



B := 



A 
W I 



and separable convex functions f : ^ g : V ^ 1 presented by comparison 
oracles, solves in time polynomial in {A, W, g{B), I, u, I, u, b, f, g) , the program 



mhi{f{Wx) + g{x) : x e Z" , Ax = 6 , i<Wx<u, I < x < u} 
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Proof. Define h : Z"+'^ ^ Z by h(x,y) := f{—y) + g{x) for all x G Z" and y E U^. 
Clearly, h is separable convex since /, g are. Now, our problem can be rewritten as 

min{Mx,Z/): (x,y) GZ"+^ \ \[^\ = [^^\^<x<u-u<y<-l}, 



and the statement follows at once by applying Lemma 2.6 to this problem. □ 



We can now provide our new theorem on generalized n-fold integer programming. 



Theorem 2.9 For every fixed integer (r, s) x t bimatrix A and integer {p, q) x t 
himatrix W, there is an algorithm that, given n, l,u E Z^, l,u E Z^"'', b E Z''"'""'^, 
and separable convex functions f : IP+^i —>■ Z, g : Z"* Z presented by comparison 
oracles, solves in time polynomial in n and {I, u, I, u, b, f, g) , the generalized problem 



mm 



n) 



x) + g{x) ■ xElT^\ = 6 , i< W^'^^x <u, l<x<u^ 



Proof. First use the algorithm of Lemma 2.7 to compute the Graver basis Q{B) of 

._ , ^^"^ 



Now use the algorithm of Lemma 2.8 to solve the problem in polynomial time. □ 



3 Multicommodity flows 

3.1 Many-commodity transshipment 

We begin with our theorem on nonlinear many-commodity transshipment. As in 
the previous section, /, g denote the maximum absolute values of the objective 
functions /, g over the feasible set. It is usually easy to determine an upper bound 
on these values from the problem data (for instance, in the special case of linear 
cost functions /, g, bounds which are polynomial in the binary-encoding length of 
the costs «£, 7e, capacities u, and demands d^, readily follow from Cramer's rule). 



Theorem 1.1 For every fixed digraph G there is an algorithm that, given I com- 
modity types, demand djj G Z for each commodity k and vertex v, edge capacities 
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Ue G Z+, and convex costs /g, : Z ^ Z presented by comparison oracles, solves in 
time polynomial in I and {d^,Ue, f,g), the many- commodity transshipment problem, 

e \ \fc=l / k=l 

I 

S.t. G Z , J2 Yl ^e=d';, J2x^,<Ue, X^ > . 

eeS+{v) e£S-{v) k=l 

Proof. Assume G has s vertices and t edges and let D be its s x t vertex-edge 
incidence matrix. Let / : Z* — Z and (7 : Z'* — Z be the separable convex functions 
defined by f{y) := Y!e=i feiVe) with ye := E1=i and ^(x) := Yfe=i EL=i dei^e)- 
Let X = (x^, . . . , x') be the vector of variables with x'' G Z* the flow of commodity 
k for each k. Then the problem can be rewritten in vector form as 



min 




+ gix) : X G Z'* , Dx'' = d^ , ^x'' <u, x>0 



fc=i 



We can now proceed in two ways. 

First way: extend the vector of variables to ') with x° G Z* 

representing an additional slack commodity. Then the capacity constraints become 
Ylk=o x^ = u and the cost function becomes f{u — xq) + (/(x^, . . . , x') which is also 
separable convex. Now let A be the (t, s) x t bimatrix with flrst block Ai := It the 
t X t identity matrix and second block A2 := D. Let d^ := Du — Yl\=i and let 
h := (m, rf*^, (i^, . . . , c?'). Then the problem becomes the (/ + l)-fold integer program 



mm 



{/ (m - x°) + ^ (x\ . . . , x') : X G Z('+i)* , yl^'^x = 6 , x > O} . 



By Theorem 2^ this program can be solved in polynomial time as claimed. 

Second way: let A be the (0, s) x t bimatrix with flrst block Ai empty and second 
block A2 := D. Let W be the (t,0) x t bimatrix with flrst block Wi := It the t xt 
identity matrix and second block W2 empty. Let b := [d^ , . . . , d'') . Then the problem 
is precisely the following /-fold integer program, 

min {/ (W^^^x) + g (x) : x G Z'* , A^'^x = b , W^^'^x <u, x > O} . 



By Theorem 



2.9 



this program can be solved in polynomial time as claimed. □ 
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3.2 Multicommodity transportation 

We proceed with our theorem on nonhnear multicommodity transportation. The 
underlying digraph is Km,n (with edges oriented from suppliers to consumers). The 
problem over any subdigraph G of Km,n reduces to that over Km,n by simply forcing 



capacity on all edges not present in G. Note that Theorem 1.1 implies that if m, n 
are fixed then the problem can be solved in polynomial time for variable number 
/ of commodities. However, we now want to allow the number n of consumers 
to vary and fix the number / of commodities instead. This seems to be a harder 
problem (with no seeming analog for non bipartite digraphs), and the formulation 
below is more delicate. Therefore it is convenient to change the labeling of the 
data a little bit as follows (see Figure [2] below). We now denote edges by pairs 



Multicommodity Transportation Problem 

Find integer I commodity transportation x of minimum f,g cost 
from m suppliers to n consumers in the bipartite digraph K^^ ^ 

Also given are supply and consumption vectors s' and d in Z', 
edge capacities Ujj and volume V|^ per unit commodity k 

consumers 




For suitable (ml,l)xml bimatrix A and suitable (0,m)xml bimatrix W 
derived from the V|^ the problem becomes the n-fold integer program 

min { f(W(")x)+g(x) : x in Z"""', A^^^x =(s', d), w(")x<u, x>o} 



Figure 2: Multicommodity Transportation Problem 
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where 1 < z < m is a supplier and 1 < j < n is a consumer. The demand 
vectors are now replaced by (nonnegative) supply and consumption vectors: each 
supplier i has a supply vector G Z'^ with s]. its supply in commodity k, and 
each consumer j has a consumption vector G Z'^ with c[ its consumption in 
commodity k. In addition, each commodity k may have its own volume Vk G Z+ per 
unit flow. A multicommodity transportation is now indexed as x = [x^, . . . , x^) with 
x^ = {x{ -^^, . . . , x-{ I, . . . , 1, . . . , where x]^ is the flow of commodity k from 

supplier i to consumer j. The capacity constraint on edge (i, j) is X]fc=i ^fc^i fe — 

and the cost is fij (eL=i ^fe^fe) + E1=i £/i,fe KJ with : Z ^ Z convex. 

As before, /, g denote the maximum absolute values of /, g over the feasible set. 



Theorem |1.3| For any fixed I commodities, m suppliers, and volumes Vk G Z_|_, there 

ij G Z+, 



is an algorithm that, given n, supplies and demands s\ G ZL, capacities Ui j G Z^ 



and convex costs fi,j,glk : Z ^ Z presented by comparison oracles, solves in time 
polynomial in n and {s'^,c\u, f , g) , the multicommodity transportation problem, 

min J2 Yl ^^'^fc + Yl Sik i^ik) 

i,j \ \ k J k=l 

I 

S.t. G Z , J2 4fc = ' Yl 4fc = 4 ' Y ^ ' 4fc ^ 

j i k=l 



Proof. Construct bimatrices A and W as follows. Let D be the (/, 0) x I bimatrix 
with flrst block Di := J/ and second block D2 empty. Let V be the (0, 1) x / 
bimatrix with flrst block Vi empty and second block V2 '■= (fi, . . . ,vi). Let A be 
the {ml, I) X ml bimatrix with first block Ai := Imi and second block A2 := D^"^\ 
Let M/^ be the (0, m) x m/ bimatrix with first block Wi empty and second block 
W2 := V^"'\ Let b be the {ml + n/)-vector 6 := {s\ . . . , s™, c\ . . . , c"). 

Let / : Z"™ Z and ^ : Z"™^ ^ Z be the separable convex functions defined by 
fiy) ■= AjiVij) with yij ■= J2[=iVkx[k and g{x) := Ei,j ELi ^l^KJ- 

Now note that A^^'^x is an {ml + nl)-vectoi, whose first ml entries are the fiows 
from each supplier of each commodity to all consumers, and whose last nl entries 
are the fiows to each consumer of each commodity from all suppliers. Therefore the 
supply and consumption equations are encoded by A^'^'^x = b. Next note that the 
nm-vector y = {yi^i, . . . , ym,i, ... , yi,n, • • • , ym,n) satisfies y = W^'^'^x. So the capacity 
constraints become W^"'^x < u and the cost function becomes f{W^^'^x) + g{x). 
Therefore, the problem is precisely the following n-fold integer program. 



min {/ +g{x) : x G Z"™^ , ^("^x = b , W^^'^x <u, x > O} 
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By Theorem |2.9| this program can be solved in polynomial time as claimed. □ 



4 Discussion 



We conclude with a short discussion of the universality for integer programming of 
the many-commodity transportation problem and the complexity of our algorithms. 

Consider the following special form of the n-fold product. For an integer s x t 
matrix D, let D'"' := A^'^^ where A is the {t, s) x t bimatrix A with first block 
Ai := It the t x t identity matrix and second block A2 := D. We consider such 



n-fold products of the 1x3 matrix I3 := [1,1,1]. Note that I3 
[3 + n) X 3n incidence matrix of the complete bipartite graph K^ n 



is precisely the 
For instance, 



[3] 



/100100100\ 
1 1 1 

1 1 1 

1 1 1 
1 1 1 

\000000111/ 

The following surprising theorem was proved in [5] building on results of [1]. (For 
further details and consequences for privacy in statistical databases see [5l El [12].) 

The Universality Theorem [5J Every rational polytope {x G : Bx = b} 
stands in polynomial time computable integer preserving bijection with some polytope 



X G 



X 



(1) 



In particular, every integer program can be lifted in polynomial time to a program 
over a matrix 1^"^''' which is completely determined by two parameters n and / only. 

Now note (see proof of Theorem 1.1) that the integer points in ([T]) are precisely 
the feasible points of some (/ — l)-commodity transshipment problem over K^ n- So 
every integer program can be lifted in polynomial time to some /-commodity pro- 
gram over some K^ n- Thus, the many-commodity transportation problem, already 
over the digraphs K^ n with fixed number 3 of suppliers, is universal for integer 
programming. So, in particular, the /-commodity transportation problem over K^^n 
is NP-hard when both n, I are variable, but polynomial time solvable for arbitrary 



fixed number n of consumers and variable number / of commodities by Theorem 1.1 
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Our algorithms involve two major tasks: the construction of the Graver basis of 
a suitable n-fold product in Lemmas |2 . 5 1 and 2.7 and the iterative use of this Graver 



basis to solve the underlying (convex) integer program in Lemmas 2.6 and 2.8 The 



polynomial time solvability of these tasks is established in [21 E]. Here we only 
briefly discuss the complexity of the first task in the special case of a digraph, which 
is relevant for the complexity of the many-commodity transshipment application. 

Let D be the s x t incidence matrix of a digraph G. Consider /-fold products 
D^'l of the special form defined above. The type of an element x = (x^, . . . in 
the Graver basis Q{D^''^) is the number of nonzero blocks x*^ G Z* of x. It turns 
out that for any digraph G there is a finite nonnegative integer g{G) which is the 
largest type of any element of any Q{D^''^) independent of /. We call this new digraph 
invariant g{G) the Graver complexity of G. The complexity of computing Q{D^'-^) 
is 0(1^^'^^) (see [2j) and hence the importance of g{G). Unfortunately, our present 
understanding of the Graver complexity of a digraph is very limited and much more 
study is required. Very little is known even for the complete bipartite digraphs K^^n 
(oriented from one side to the other): while glK^^s) = 9, already glK^^^) is unknown. 
See [1] for more details and a lower bound on g{K^^n) which is exponential in n. 
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